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Abstract 

The concept of Universal Shower Profile is used to characterize the average behavior of high en- 
ergy cosmic rays. The shape variables contain important information about composition. They are 
independent of the primary cross-section by construction, but affected by other hadronic parame- 
ters, like multiplicity. The two variables give access to the average nuclear mass of the sample and 
their compatibility serves as a test of hadronic models. 
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1. Introduction 

The development of high energy cosmic ray showers in the atmosphere gives indirect informa- 
tion on the primary particle. It is, however, governed by very high energy hadronic interactions 
for which there is no unique description, only different phenomenological models that try to ex- 
trapolate the accelerator data obtained at lower energies and different phase space regions. 

In general the maximum number of particles and the depth at which that maximum occurs 
are used to characterize the showers. The first is a measurement of energy, while the the second 
statistically distinguishes primary particles and depends on the cross-section of the first interaction. 
More information is however hidden in the details of the shower longitudinal profile shape [QH^S] • 

Electromagnetic (and muonic fl) longitudinal shower profiles are well described by the Gaisser- 
Hillas function, which can be parametrized in N' = N/N max and X' = X - X max as: 



, . RX'\ R 2 / X' . 
N = 1 + -H ex P -777 (D 



With the above parametrization, introduced in fij^, the profile can be recognized as a Gaussian, 
with width L, with an asymmetry introduced by non-zero values of R. Notice that after translation 
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to the maximum, the variations due to the point of first interaction disappear 

We will start by introducing the Universal Shower Profile (USP), showing that the shape of 
the energy deposit longitudinal profile is almost universal. Both shape variables, L and R, tend 
to different values according to energy, primary particle and hadronic interaction model. We will 
then study the construction of the USP for a given energy, regardless of the possible mixed compo- 
sitions. In fact, both parameters are a function of the average nuclear mass, with different energy 
evolutions. We will find that the different hadronic interaction models define specific possible re- 
gions for the two parameters, which can then be used to exclude or constrain present models or 
help in the construction of new ones. 



2. Shower Profile Universality 

The universality of the shower longitudinal profiles, when expressed in N' = N/N max and 
X' = X- X max was discussed in yj], where it was shown that L and R could give extra information 
on the shower properties. Here we concentrate on what can be learned from the average showers. 
In fact, the stability of the above variables within showers initiated by similar primaries means 
that by collecting only a limited number of events, important information can already be obtained, 
with a significant reduction of statistical uncertainties. 

The average shower profile is not exactly described by a Gaisser-Hillas and so the values of 
L and R can vary slightly according to the X' interval considered in their determination. In order 
to increase the available data in a real experiment, the total X' interval should be minimized. 
However, for a good determination of the width, L, a minimum range around the maximum is 
needed, while for R the rising tail needs to be includedyj]. From now on, we will always consider 
the range of X' e [-350,200] g/cm 2 . 

Figure Q] shows the build up of the USP variables, with the collection of events. The energy 
deposit shower profiles were generated with CONEXlfiiS], using the QGSJet-II.03[Q, 0] model, 
for different primaries at fixed energy of log(£'/eV) = 19. After around one hundred events, the 
average shape becomes stable, with a good separation for the different primaries, in terms of both 
variables. 

The USP variables evolve with energy, and figure [2] shows the results for proton and iron 
primaries of different energies. Average values of L and R are reached when different energies 
are mixed. In general, data analysis will mix a range of energies, and care must be taken in the 
interpretation of the results. However, it can be seen that even for a wide interval of log(£7eV) e 
[18.5, 19.0], light and heavy primaries can still be distinguished. 



3. USP and Mass Composition 

The values of L and R have an almost linear dependence with the logarithm of the nuclear 
mass, ln(A), as can be seen in fig. [3j There, in addition to the pure samples, also different mixtures 



Although changes in the first interaction point affect the atmospheric profile crossed by the air shower, the cor- 
responding effect on the profile shape in gem 2 is negligible when compared to the differences between the several 
hadronic interaction models used for the shower development, for example. 
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Figure 1: Shape parameters L (left) and R (right) as a function of the number of events used to get the average 
shower. The showers were generated using QGSJet-II as the high energy hadronic interaction model. Color codes 
show different primaries atlog(£'/eV) = 19.0. 
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Figure 2: Shape parameters L (left) and R (right) as a function of the number of events used to get the average shower. 
The showers were generated using QGSJet-II as the high energy hadronic interaction model. Darker colors show fits 
to the USP of proton and iron showers at both log(£'/eV) = 18.5 and log(£'/eV) = 19.0; light colors correspond to 
the USP obtained by mixing events generated at two energies evenly. 
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Figure 3: Shape parameters L (left) and R (right) as a function of (ln(A)}. The showers were generated using QGSJet- 
II as the high energy hadronic interaction model. Color codes show different combinations of primaries, generated at 
log(£/eV) =19. 



of primaries are presented. Even in these cases, the stability around an average shape is easily 
achieved. 

The values of both parameters give two measurements of (ln(A)), with small variation accord- 
ing to the mixture of the sample. A pure sample of helium primaries will have an only slightly 
lower value of R and higher value of L than a mixture of proton and iron primaries with the 
same average logarithmic mass. Since the values of both parameters correlate with (ln(A)), a well 
defined region can be expected in the L, R plane, as shown in fig. |4l 

Notice also that the same analysis can be extended to include X max , bearing in mind that in that 
case there is an influence of the primary cross-section: a big change of the cross-section (as, for 
example, suggested in ref. DSD) would then be apparent as a consistent reading of (ln(A)) between 
L and R, with an inconsistency for (X max ). 



4. Primary composition and hadronic parameters 

The different values of the shape variables L and R according to primary composition, energy 
and hadronic interaction models can be used to understand the hadronic parameters governing 
them. While the new observables are insensitive to the point of first interaction, they can be deeply 
changed by the most usual processes in cascade development. 

The correlated evolution of L and R is characteristic of the hadronic interaction model. Figure 



|5]shows the comparison of four different high energy hadronic interaction models: QGSJet01c[10], 



SIBYLL2.1H1 10 and EPOS1.99UJ], in addition to QGSJet-II.03, for showers at log(£/eV) = 19. 
The measurement of one point in this plane can be used to check the compatibility of each model 
with experimental data. 

It can be seen, for example, that EPOS has larger L values for proton showers than the others; 
it is even more interesting to note that QGSJet-II spans a range of the L, R plane which is quite 
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Figure 4: Relation of the shape parameters of the universal shower profile with the (ln(A)}. The showers were 
generated using QGSJet-II as the high energy hadronic interaction model. Color codes show the (ln(A)) for different 
combinations of primaries, generated at log(£'/eV) = 19, and crosses mark the results for pure composition samples. 



different from the other three models, namely in what concerns the larger R values for heavy 
nuclei. This might reflect the fact that multiplicity is much higher in this model. High values of 
R characterize heavy primaries, which also have larger multiplicities in the first interactions, and 
the same kind of physical characteristics that enable the distinction of primaries is expected to be 
present between hadronic interaction models. 

Figure [6] shows the evolution of the parameters as a function of energy, for the different models. 
Notice that the electromagnetic energy is given by E em = ^FhiL ^|| and L should be mainly 
determined by the primary energy. There is a relative change of the order of 5% between primaries, 
and of the order of 1% between models. On the other hand, R migth be related to multiplicity, 
with large differences between primaries. We highlight these trends by using the scaled energy, 
logCE/A), inspired in the superposition model for showers initiated by nuclei. The scaling is 
reflected within each model, with all primaries and energies lying approximately in two lines: one 
for QGSJet-II and another one grouping all other models. 
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Figure 5: Shape parameters (L, R) of the universal shower profile , generated at log(£'/eV) = 19, for the different 
mass composition combinations in several models QGSJet-II.03, SIBYLL2.1, EPOS 1.99, QGSJetOlc. 
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Figure 6: Shape parameters L (left) and R (right) as a function of energy for four different primaries (p, He, N and Fe) 
at log(£/eV) = 19.0 in all models, shown in colored dots, and also log(£/eV) = 18.5 in QGSJet-II and QGSJetOlc, 
shown as squares and log(£'/eV) = 18.5 in QGSJet-II, as triangles. In the case of R, the energy is scaled by mass 
number, to highlight the hidden dependency (see text). 



5. Summary and Conclusions 

A Universal Shower Profile can be constructed with a limited number of events to characterize 
the average behavior of high energy cosmic rays. By constructing an average longitudinal shower 
profile shape, the statistical errors associated to the detection of single events are significantly 
reduced. The average shape can be described with the variables L and R, which contain important 
information on composition. 

While being independent of the first interaction point by construction, the Universal Shower 
Profile shape keeps information about the other parameters governing shower development which 
are different among primaries and models, namely multiplicity. Within each model, L and R 
give two measurements of (A), with a reduced dependence on the primary cross-section. The 
consistency between both is a powerful test of the modeling of the hadronic interactions. 
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